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Abstract 

This is the third paper of a series devoted to a systematic study of QCD correlation 
functions in the framework of an instanton model for the QCD vacuum. In this paper we 
concentrate on the quark-quark (diquark) and three quark (baryon) channels. We have 
found that the quark-quark interaction resembles the one between quarks and antiquarks. 
Similar to the pion and rho channels, the interaction in the scalar isospin 1 = and vector 
1 = 1 diquark channels is completely different: the former has a much stronger attractive 
interaction. As a consequence, the SU{3) octet (nucleon) and decuplet (delta) correlators 
are also found to be qualitatively different. Using a complete set of all available correlation 
functions, we determine masses and coupling constants for the nucleon and delta. Our 
results agree surprisingly well with the first lattice data on point to point correlators. 
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1. Introduction 



This is the third paper of a series devoted to a detailed investigation of the QCD 
point-to-point correlation functions^ in the framework of the Random Instanton Liquid 
Model (below RILM). After considering the propagation of quarks and quark-antiquark 
(mesonic) states we now consider the baryon sector. 

This model of the QCD vacuum was proposed by one of us @] in 1982. Using a variety 
of arguments, it was suggested that the density of 'tunneling events' (instantons and 
antiinstantons combined) is about uq = 1 /m~^, while their typical size is po = 1/3 /m. 
These two numbers were shown to reproduce two global properties of the QCD vacuum, 
the gluon and quark condensates. 

The investigation of correlation functions in the RILM was started in reference for 
the pseudoscalar channels. Not only does the model explain why the OPE-based sum 
rules fail in the case of pseudoscalars but it was also shown to reproduce deviations 

from asymptotic freedom at small distances in a quantitative way. 

During the past decade significant efforts were devoted to derive the instanton picture 
of the QCD vacuum from first principles. First variational [§] and later numerical 0, 
p!0| methods were applied to this problem, with the hope to develop a self-consistent 
theory of interacting instantons. In spite of these efforts, a clear cut derivation from first 
principles is still missing. Recently, we have performed further studies along these lines, 
substituting the rather arbitrary trial functions used previously by the 'stream-line' set of 
configurations found in |Tl|. These studies have shown that we still lack an understanding 
of how exactly one should treat 'a bottom of the valley', the configurations with a very 
close instanton-anti-instanton pair. We are now preparing a separate publication, in which 
the current understanding of this problem will be discussed. 

In the RILM we simply assume that except for the size, which we keep fixed, the 
distribution of the collective coordinates associated with the instantons and antiinstantons 
is completely random. We consider the RILM as the simplest possible model of its class 
detailed review on QCD correlation functions recently appeared in 
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and a benchmark for any more involved calculation with 'correlated instantons'. Before 
studying more complicated ensembles, it is important to understand the predictions of 
the random model for physical quantities. Our previous studies of the propagation of 
a single quark and quark- ant iquark pairs (mesonic correlators) in such a vacuum 
have significantly extended those reported in [|12[: the box volume was increased by a 
factor 20, the correlation functions were traced down by a few orders of magnitude, etc.. 
Furthermore, the results were shown to be in quantitative agreement with experimental 
data, which was not expected for this simple model. Experimentally known correlators 
(and eventually, masses and coupling constants of the mesons) are reproduced, in some 
cases literally inside the uncertainties of the calculation! What is equally exciting is the 
fact that we also observe good agreement with the results of recent lattice calculations 



T3[| , in spite of the quite different approaches and approximations used. 



As this paper is the first exploratory study of baryon correlators in the instanton 
vacuum, the use of the random model appears to be well justified. Its primary goal is to 
check whether this model leads at least to some baryonic bound states, in the form of a 
resonances well separated from a multi-particle continuum. 

Before we start, let us make a few general remarks about the baryonic channels. 
It is not surprising that quark-quark forces are even less understood than the quark- 
antiquark ones: ordinary (made of light quarks) baryons are not simple two-body systems. 
The situation is somewhat simpler if one of the quarks is 'sterile' (very heavy) |T^, but 
unfortunately the experimental data on heavy-light baryons are still rather incomplete. 



A direct manifestation of the 't Hooft instanton-induced interaction |T5| (which pro- 
duces so spectacular effects for the scalar and pseudoscalar mesonic correlators) is also 
expected in the quark-quark interaction. In this case, it appears as a strong attraction 
in the scalar diquark channel [jT^. As suggested in [T^, this interaction can naturally 



explain all spin splittings in the octet and the decuplet, describing the spectra at least as 
good as traditional explanations using gluomagnetic forces between 'constituent' quarks. 
Thus, our second major question is whether baryons belonging to the two lowest SU{3) 



multiplets, the octet and the decuplet, are similar (their differences being due to relatively 
small perturbative effects), or whether they are completely different. In this context we 
not only want to consider the N — A mass splitting, but also make a comparison of the 
coupling constants and the general shape of the correlation functions. 

The paper is structured as follows: in section 2 we introduce 'diquark' correlators, 
which are compared to the corresponding mesonic channels and then used (in section 3) 
to derive properties of heavy-light baryons. In sections 4 and 5 we deal with the nucleon 
and delta correlators, as representatives of the SU{3) octet and decuplet baryons. A 
comparison with QCD sum rules and lattice results, as well as a general discussion is 
contained in sections 6 and 7. Conclusions are summarized in section 8. 



2. Diquarks 

In this section we study correlation functions of two-quark currents. As suggested in 
[1^ , such currents can be thought of as part of a baryon current in which one of the 



quarks is 'sterile', or very heavy, so that it does not move at all. We will come back to 
this interpretation in the next section. 

The generic diquark current is defined by 

J? = e'^'^q'^iCTUrql,, (2.1) 

where C is the charge conjugation matrix and F is one of the possible gamma matrix 
structures. Furthermore, denotes an isodoublet quark spinor with color index a and 
spinor index a. The isospin matrix r is given by r = r2 for an isosinglet diquark and 
by r = T2T for an isovector diquark. Note that the isospin wavefunction is determined 
by the symmetry of the Dirac matrix CT and the overall antisymmetry with respect 
to exchanging the two quarks. In particular, F = {1,75,7^75} corresponds to isoscalar 
diquarks whereas F = {7^, a^^} gives isovector diquarks. 

Note that the parity of the Dirac matrix CF is opposite to that of the matrix F: 
the internal parity of a qq pair is positive while the one of the qq system is negative. 



For example, if F = 75, the diquark is a scalar. This well known fact makes diquark 
notations somewhat ambiguous, one can either (i) specify the real quantum numbers of 
the di quarks, or (ii) indicate which gamma matrix was used in the current. We use the 
latter notation (e.g. call the F = 75 case 'the pseudoscalar correlator'), since in this 
notation the comparison between diquark and meson correlators is more straightforward. 
We hope that this will not confuse the reader. 
We study diquark correlators defined by 

nr(x - y) =< j^{x) ij^iy))^ >, (2.2) 

where the hermitian conjugate current {j^V is given by 

(Jr)^ = e^'^ti^CUrr,. (2.3) 

The correlator can be expressed in terms of the quark propagator resulting in 

nr(x -y) = TT [S''%x - y)rCS^''\x - y)Cr] - Tr [s^%x - y)rCS^^\x - y)Cr] ,(2.4) 

where the transpose is with respect to the Dirac indices. At short distances we can 
evaluate the correlator using free quark propagators. Deviations from the free quark 
propagation have been studied using the operator product expansion (OPE), further sim- 



plified by the so called 'vacuum dominance' approximation For baryonic currents it 
was implemented |]19[ in a very simple way: the propagator is given by the sum of a free 
propagator and a quark condensate. As discussed in more detail in [|I^], this approxima- 
tion corresponds to a 'distance dependent constituent quark mass' ttiq = (7r^/3)r^| (gg) | 
which is generated as the quark propagates in the chirally asymmetric QCD vacuum. The 
corresponding propagator reads^ 

^S{r) = -^70 -^<qq>. (2.5) 

As we have shown in the case of the mesonic correlators, this approximation is quite useful 
in understanding the experimental trend of the various correlation functions, although it 



In what follows we will always work in euclidean space-time. Physically there is no preferred direction, 
so we can always choose x — y to point in the euclidean time direction. 



is far from being quantitatively correct. For the diquark correlators introduced above we 
find 
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(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 
(2.6e) 

where F = {S, P, V, A, T} labels the gamma matrices in the diquark channel and np 
denotes the free correlation function. Note that these estimates for the short- distance 
behaviour of the diquark correlation functions exactly coincide with the corresponding 
approximation for the mesonic correlator in the same gamma matrix channel. This is a 
consequence of the fact that in the vacuum dominance approximation the propagator is 
diagonal in color space, which is not correct in general. 

We now come to a physical interpretation of the diquark correlators. Since diquarks 
are color antitriplets, they cannot appear as asymptotic states. Their correlation functions 
are gauge-dependent and unphysical, unless either (i) a gauge fixing is implied or (ii) the 
additional color matrix U = Pexp{ig dx^A^j^t"- /2) is included. If the path C is a segment 
of a straight line going in the time direction, it corresponds to a static third quark: thus 
we naturally come to light-heavy baryons, to be discussed in the next section. 

The specific feature of instanton models is that the matrix U ^ 1, because heavy 
quarks interact very little with instantons. We will therefore neglect the effects of U in 



our following analysis. More quantitatively, it was shown in that the gauge string U 
produces a potential of the order of 50 MeV, which is well inside the accuracy of the present 
calculations. Thus, inside instanton-based vacuum models (which ignore confinement etc.), 
the diquark correlators defined above are 'practically gauge invariant '0. In such a spirit, 
^ The reader who does not Uke such an approximation, may skip our quaUtative discussion of di- 
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one can define the coupling of the currents introduced above to 'diquark states' 



<mp\^s,P> = ^^'^71^"'"' (2.7a) 

Sjab 

<0\fv,AMv,A> = 3v,A^e,e-'''^, (2.7b) 

sab 

<0\jtr,AMs,P> = fv,A^Q,e-"^'', (2.7c) 

where |0p > denotes a diquark state with momentum and color index a that has 
the quantum numbers of the corresponding current and is the polarization of a vector 
diquark. Proceeding this way we effectively treat diquarks just like mesons. In partic- 
ular, we also allow for mixing of vector and scalar diquarks in the vector channel. This 
phenomenon is analogous to the well known tt — oi mixing which is observed in the axial 
vector channel. 

In the following we want to describe the diquark correlation functions in terms of 
the physical intermediate states which can contribute. For this reason we represent the 
euclidean Fourier transform of the correlation function by a dispersion integral 



/ 



d'xe'"nrW = ir<i. '°'"'-W . (2.8) 

n Jo q'^ + s 



where pr(s) = ;^Imnr(s) is the spectral density associated with the corresponding corre- 
lation function. The spectral density can be expressed as 

pr{s = e) = {2T,fY.^\k - qn) < 0|jr(0)|n >< n|j^(0)|0 > , (2.9) 

n 

where \n > is a physical intermediate state with momentum 

Similar to the procedure in the meson sector we approximate the spectral function by 
the contribution of one or two diquark resonances and a free quark continuum starting 
at a threshold Sq. This means that for long distances the correlator is dominated by the 
contribution of the lowest physical state whereas for short distances it can be described 
in terms of the propagation of free quarks. In the scalar and pseudoscalar channels we 



quarks and proceed directly to our discussion of nucleon and delta correlators, which are explicitly gauge 
invariant. 



get 

3 

Ps,p{s) = gl,pS{s - mlp) + -^sQ{s - So), (2.10) 

where the first term is the resonance contribution and the second term represents the 
continuum. On purely dimensional grounds we have ps,p{s) ~ s for large invariant masses. 
The corresponding coefficient can be calculated from the discontinuity of the free quark 
loop diagram appearing in eq.( |2.4| ). 

Inserting the model eq. ( P^.10[ ) for the spectral function in the dispersion relation ( |2.8| ) 
and transforming back to coordinate space the diquark correlator is given by 

3 

n5,p(r) = glpDims,p,T) + — / dssDiV^s^r) , (2.11) 

where D(m, r) = m/ (47r^r)il'i(mr) is the euclidean space propagator for a massive scalar 
particle. Proceeding in the same way in the vector channels, we find 

3 

ny,A(r) = 3gy AD{mv,A, r) - mlpfy .D{ms,p, t) + — ds sD{^/s, r) , (2.12) 

ZTC Jso 

In the tensor channel, finally, there is no continuum contribution. The dirac tensor a^^, 
has components O"oi which transform like a three-vector and eij^cy-jk which transform like 
a pseudovector. This means that vector as well as axial vector diquark states can con- 
tribute. Since the contributions from these two states are hard to distinguish in practice 
we describe the tensor diquark correlation function using a single diquark resonance with 
mass rriT and coupling qt 

IiT{r)=glD{mT.T). (2.13) 

We have extracted the diquark masses and coupling constants by fitting the parametriza- 
tions eq.( p.ll[ - p.l3| ) to the diquark correlation functions measured in the random instanton 
liquid. The correlation functions together with the fits are shown in figure 1 and the results 
of the fit are summarized in table 1. 

Looking at the various diquark correlation functions one immediately realizes that 
there are very pronounced qualitative differences between the different channels. In par- 
ticular, the scalar case (F = 1) has very repulsive interaction and shows no indications for 
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a 'diquark resonance' whereas the pseudoscalar one (T = 75) shows a huge enhancement^, 
or strong attractive interaction. Determining the parameters of this resonance we find 



mp = 420 ± 30 MeV 



gp = 0.23 ± O.OlGeV^. 



(2.14) 



The dashed curve in figure 1 shows the result of the 'vacuum dominance' approxi- 
mation, see 0. As for the mesonic channels, it reproduces the qualitative differences 
between the scalar and pseudoscalar channels but fails to give a quantitative description 
a distances larger than 0.3 fm. 

In the (r = 7^) vector diquark channel the correlation function remains close to 
the perturbative one out to fairly large distances. Again, the situation is similar to the 
superduality phenomenon observed in the rho meson channel. The resonance parameters 
in the diquark channel are 



In the (r = 7^75) axialvector diquark channel, we can have mixing with the light pseu- 
doscalar diquark. Since the contribution from this state alone can already explain the 
measured correlation function, we have no evidence for an axialvector 'diquark state' in 
the RILM. In the tensor (F = cr^iy) diquark channel we have fitted the correlator with a 
single resonance at mx = 570 ± 20 MeV, which is intermediate between the light pseu- 
doscalar and heavy vector diquark resonances. 

Of great interest for hadronic models which invoke diquark clustering in order to ex- 
plain the observed spin splitting between the octet and the decuplet is the mass difference 
my — mp c::^ 520 MeV. In a simple quark-diquark model p!7| , pO|] we have — mj^ = 
I (my — mp). Comparing with the experimental splitting mA — m^ = 293 MeV we do in 
fact get a very good agreement. We will further comment on the quark-diquark picture 
of the baryons in our discussion of the nucleon and delta correlation functions. 

■^Which is still not as strong as for the pseudoscalar mesons, of course, which have to be massless 
due to Goldstone theorem. However, it is quite comparable with the effect observed for scalar isoscalar 
'sigma' meson, which also has a comparable mass. 



mv = 940 ± 20 MeV, 



gv = 0.24±0.01GeVl 



(2.15) 
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3. Heavy-light baryons 



Adding a heavy quark Q to the diquark currents considered in the last section we can 
construct currents carrying the quantum numbers of heavy baryons. Correlation functions 



of these currents were first considered in |14| in the framework of the QCD sum rules 



Recently, further analysis along these lines was performed in [21, where somewhat 
different (non-relativistic) currents have been used. 

The heavy baryon currents have the general structure Jsga = jr(J"Q)a where jr is a 
diquark current, Q is a heavy quark spinor with dirac index a and T' is a dirac matrix. 
Out of the five possible currents, three couple to the Aq while the two remaining ones 
have the quantum numbers of the Eg: 

Jaq = e-'%u''C-f,-f5d')YQ', (3.1a) 
Ji^ = e"^^(M«C7^u'')757^Q' = e-'%u-Ca^,u'')^,a^'''Q'^ . (3.1b) 

In the following we will consider correlation functions constructed from these currents 

^ap{x -y) = i < JBQa{x)JBQl3{y) > ■ (3.2) 

It is instructive to consider the limit in which the mass of the heavy quark is taken to 
infinity. In this case the correlation function simplifies to 

^afi{x -y) = {T'iSQ{x - y)T')af3 ^r{x - y), (3.3) 

where V' = 7or'^7o and 

^S<,{x -y) = l±^e(r)5^(f - ,1 (^) ' e""^ (3.4) 

denotes the heavy quark propagator in the nonrelativistic (static) limit. Using this ap- 
proximation the spin structure of the correlation function becomes trivial and all the 
dynamical information is contained in the diquark correlator IIt{x — y). 



As in the last section we want to describe the Euchdean correlation function (|3.2| ) in 
terms of physical intermediate states. This is achieved by means of the spectral decom- 
position of the correlator 

roo 

^apir) = dspa(3{s)D{y/s,T) (3.5) 

JO 

Again we approximate the spectral function by the sum of a resonance and a continuum 
contribution, which in this case results in 

s 

+ - ^o)'/'^(s - so)((7 ■ q)ap + Mq5^p). 

The coupling constant fs is defined by 

< 0| Jb^„|5q >= fB^Bi^\ (3.7) 

where Bj^^ is a Dirac spinor with spin s. The numerical factor multiplying the continuum 
contribution follows from asymptotic freedom. 

It is convenient to measure the energies with respect to the heavy quark mass Mq. 
When we write s = {Mq + cj)^, the non-relativistic heavy quark propagator factorizes 
from the spectral function, and the Euclidean diquark correlator is given by 

R 1 r°° 

iir(rj = — ^exp(-ei3Qr) + -— / du exp{-ujT) , (3.8) 

where e Bq is the heavy baryon energy (minus the rest mass of the heavy quark) . As above 
we have introduced a continuum integral starting at a threshold energy ujq- 

The heavy baryon parameters extracted from this non-relativistic interpretation of 
the diquark correlation functions are given in table 2, together with the predictions from 
QCD sum rules. 

The Aq results were obtained from the pseudoscalar diquark channel, since this current 
has a simple non-relativistic limit. We find 

€aq = 760 ± 30 MeV /aq = 0.052 ± 0.005 GeV^ . (3.9) 
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The parameters of the Eg baryon can be determined from the F = 7^ or F = a^j^p diquark 
channels. Using the F = 7^ diquark channel we get 

e^^ = 892 ± 30 MeV, /^^ = 0.011 ± 0.002 GeV^ (3.10) 

while the fit to the tensor channel gives a somewhat larger massQ 

el = 1006 ± 30 MeV, /J = 0.044 ± 0.008 GeV^ . (3.11) 

Of particular interest is the Eg — Aq mass splitting which we find to be J^j-Eq ~ ""^Aq = 
135 MeV or 246 MeV, depending on the current used in the fit. This number should 
be compared with m^^ — = 170 MeV for charmed baryons. Unfortunately, the 
corresponding value for bottom baryons is not known experimentally. 

It is instructive to compare these results with those that have been obtained from 



QCD sum rules. In the original paper |]14| the predictions areQ = 700 MeV and 
/aq = 0.02 GeV'^ together with a mass splitting tt^Sq — "i^Aq = 400 ± 250 MeV. On the 
other hand, the values obtained in |2| are eAg = 780 MeV and Jaq = 0.022 ± 0.005 GeV^ 
which is very similar, but there is no Sg — Ag splitting! Let us make a few remarks 
concerning these results: 

(i) The two OPE-based papers give similar predictions for the mass and the coupling 
constant of the Ag when the pseudoscalar current is used. While the predicted mass 
agrees with our fit, the coupling constant does not. This is not surprising: the large 
instanton-induced attractive interaction which leads to the strong enhancement in the 75 
channel is not reproduced by the OPE. 

(ii) Extracting the Eg parameters and the Ag — Sg splitting from the OPE is very 
subtle: in fact the 'vacuum dominance' approximation (inclusion of the quark condensate 
term in the propagator) does not lead to an unambiguous result. 

Indeed, compare for example the 7^ and 7/^75 correlators in order to determine the 

Sg — Ag splitting. In this case the corrections up to order (qq)"^ read 1 ± (gg)^r^7r^/18, 

^ Naturally, the two coupling constants , have different definitions and should not be equal. 
^There was a misprint in the paper: The value of /aq is one order of magnitude larger than indicated. 
The published version corresponds to the r.h.s. of the sum rules 100 times smaller than the l.h.s. 
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respectively, so that one may think there is a big sphtting. However, if one considers the 
non-relativistic hmit of the currents (see ref. [0), using the 7 di quark for the nonrelativis- 
tic Eg current and 7075 for the Aq, one finds identical OPE expressions 1 + (gg)^r^7r^/36 
in both cases, leading to zercJ^ mass splitting! 

(iii) Our fits for the Eg mass in the V and T channels are different, but not very much 
so. A Sq — Aq splitting on the order of 200 MeV is definitely predicted. What is maybe 
even more important is that we see a qualitative difference between the Eg and the Aq : 
both V and T correlators are comparable to the perturbative ones at distances about 1 
fm, while the correlator in the 75 channel is several times larger. 



4. Nucleon correlation functions 

After studying the quark-quark interaction we now proceed to correlation functions 
of baryonic (three quark) currents. The Euclidean correlator of two spin 1/2 nucleon 
currents is defined as 

T^ap{x -y)=i< J^{x)J^{y) >, (4.1) 

where are the spinor indices of the nucleon currents. Using Lorentz and parity 
invariance one can show that the correlator can be decomposed in terms of only two 
independent Dirac structures 

na/3(^ -y)= ni((x - yf){^ ■ {x - y))^p + U^iix - yf)5^p. (4.2) 

Again we will consider the correlation function in the euclidean time direction. Only two 
independent correlators remain for a given current, namely Tr(n(r)7o) and Tr(n(r)). 

Due to isospin symmetry it is sufficient to consider only one charge state of the nucleon. 
A proton current is constructed by coupling a rf-quark to a uu-diquark. The diquark has 
the structure eabcUbCTuc which requires that the matrix CT is symmetric. This condition 
^The contribution of the gluon condensate also turns out to be identical. 
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is satisfied for tlie V and T gamma matrix structures. Only two possible currents that 
have positive parity and spin 1/2 can be constructed [jl9|. These so called loffe currents 
are given by 



Vi = ^abc{u''Cj^u'')j5-ff,d'', (4.3a) 

In total, there exist six different nucleon correlators: the diagonal r/iz/i, and off- 
diagonal r]ifj2 correlators, each contracted with either the identity or 79. They are defined 
as 



nf = 


^ < Tr(r7ir/i) >, 


(4.4a) 


n^ - 


i < Tr(7or/ir7i) >, 


(4.4b) 


ttN _ 
113 — 


i < Tr (7/2^72) >, 


(4.4c) 


nf = 


1 

- < Tr(7or/2r72) >, 


(4.4d) 


nf = 


^ < Tr(r7ir/2) >, 


(4.4e) 


nf = 


^ < Tr(7or7ir/2) >, 


(4.4f) 



Again we can discuss the qualitative features of these correlation functions using the 
vacuum dominance model (|2.5|) . The result is 



n^ Vr^ TT^ 

W = T^K^"^)I^' + ^K5-5)|V, (4.5a) 
n^ ir'^ 

n|o = ^ + ^\m'r', (4.5b) 

nf 

nfo 
n^ 

^ = Y^mir', (4.5e) 

§0 = (4.5f) 



1, (4.5d) 
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In order to compare these functions with experimental information we have to analyze 
the correlation functions in terms of physical intermediate states. The nucleon coupling 
to the two loffe currents is given by 

<0|r/i,2(0)|iV(g,.)> = X^^,^l^^uiq,s)e-^^^, (4.6) 

where \N{q,s) > denotes a nucleon state with momentum q and spin s, u{q,s) is the 
corresponding free fermion spinor (normalized as in [^) and is a coupling constant. 
The first two correlation functions introduced above are related to the first loffe current. 
Using the simple nucleon resonance plus continuum model introduced in section 2 we 
have 

/ \ 2 1 /""^ 

nf(r) = {X^)mND{mN,T) + —\{qq)\J^^dssD{^/^,T), (4.7a) 

n^(r) = -(Af) D'{m^,T)-^J^^ dss'D\^s,T), (4.7b) 

where D{m,T) is the euclidean space propagator of a scalar particle of mass m and 
D'{m,T) = d/{dT)D{'m,T). For the second correlation function 11^ the continuum corre- 
sponds to the propagation of three free quarks and can be calculated from the disconti- 
nuity of the quark loop diagram contributing to the correlation function of the two loffe 
currents. 

In the case of Hi there is no such contribution. In order to improve the description at 
short distances we have added the discontinuity of the diagram in which one of the quarks 
interacts with the quark condensate. In the language of the operator product expansion 
this means that we have calculated the spectral function from the imaginary parts of 
the diagrams which determine the coefficients of the unit and (gg)-operators. In a more 
physical picture one can check that, at least in an average sense, the continuum part of 
Hi represents the contribution of the vrA-continuum and higher nucleon resonances. 

The correlation functions 11^4 are related to the second loffe current. Proceeding as 
above they can be written as 

nf(r) = {X^ymND{m^,T), (4.8a) 
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(4.8b) 



The last two nucleon functions are related to off-diagonal correlators between the two 
loffe currents. They are given by 



The measured nucleon correlation functions together with the physical parametrizations 
introduced above are shown in figure 3. The dashed lines show the result of the vacuum 
dominance approximation. Nucleon parameters determined from fitting the correlation 
functions in the RILM are given in table 3. 

The most important feature of the nucleon correlation functions is the strong en- 
hancement over the free correlator which can be seen in and Uf. This behavior is 
reminiscent of the pseudoscalar diquark channel which shows a similar strong attractive 
interaction. Somewhat surprisingly, for the nucleon case this feature is semi-quantitatively 
reproduced by the vacuum dominance approximation. 

This enhancement in the calculated correlators is nicely described by a nucleon con- 
tribution to the correlation function. Fitting the nucleon parameters we get 



niN = 960 ± 30 MeV, Af = 0.031 ± 0.001 GeV^ = 0.080 ± 0.004 GeV^ (4.10) 



with a surprisingly accurate mass value. What is even more important, the simple 'nucleon 
pole plus continuum' model gives a very good simultaneous description for the complete 
set of correlation functions. 

This agreement is particularly good for correlators involving the first loffe current, 
while it is somewhat worse for the other ones. Note that the correlation functions and 
(which show deviations from the simple 'nucleon pole plus continuum' model) involve 
at least two quarks which have to flip their chirality. It might be that the RILM has 
problems reproducing those amplitudes. An alternative explanation (suggested by loffe 
in his original QCD sum rule analysis) might be that the second loffe current couples 
more strongly to some excited states of the nucleon. 




(4.9a) 




(4.9b) 
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It is important to note that the simple random instanton model, even without con- 
finement, does create a nucleon hound state. The presence of a bound state is indicated 
by the fact that we can describe the correlation functions in terms of a single pole, clearly 
separated from the free quark continuum starting at i?o — 2 GeV. 

In order to strengthen that point we have also considered other models for the nucleon 
correlation functions. The simplest one is a 'constituent quark model'. One may try to de- 
scribe the correlator in terms of the propagation of three uncorrelated 'constituent' quarks. 
For their masses we take the coordinate dependent effective mass mq = (vr^/3)r^| (gg) p 
at short distances and mq = 300 MeV at large distances. The result is shown in figure 
5 for the case of 11^. The conclusion is obvious: the constituent model clearly can not 
explain the shape and magnitude of the correlation function. 

A somewhat less trivial model for the correlator is the quark-diquark description. For 
this purpose we have Fierz rearranged r/i in such a way as to make the scalar and vector 
diquark content explicit |jT9|. We can now use the measured diquark correlators which 
show a strong enhancement in the scalar diquark channel, and describe the propagation 
of the third quark using the constituent quark propagator introduced above. Again, the 
result is presented in figure 5. The corresponding curve is much closer to the data than 
the pure constituent model, but it still fails to reproduce our numerical results. 

We conclude form these comparisons that the RILM not only produces substantial 
attraction in the scalar diquark channel, but also describes the nucleon as a true three 
quark bound state. While the first feature is essentially a consequence of the lowest order 
t'Hooft interaction, the presence of a nucleon pole is a very non trivial consequence of the 
RILM. 

5. Delta correlation functions 

In the case of the A— current, isospin symmetry allows us to restrict ourselves to 
the A"*"^. In this case all quarks have the same fiavors and with the help of a Fierz 
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transformation one can show [jl9[ that there exists only one independent current which is 
given by 

= eaUu'^CjyX. (5.1) 
However, the spin structure of the correlator 

nt;a/3(^-l/) = t<jtix)J^^{y)> (5.2) 

is much richer. As before the correlator is related to the spectral function by 

nj^;a/3(a;) = / dsPf,^,a(^D{^/s,x). (5.3) 



For fixed values of fi and u we again have a two-spinor spectral function which, as seen 
in section 4, reduces to two invariant structures which can be projected out by either 
contracting with the identity or 79. The /xi/— indices transform as a Lorentz vector. Five 
different vector structures are possible 

(5.4) 

which adds up to ten different possibilities for the correlator. However, the delta current 
satisfies the Rarita Schwinger constraint 7^ J)f = which allows us to express the spectral 
function p^u;ai3 using only four independent form factors 

+ a2{q^) ilf^luq - 16^^g + 27^g^ + 47^g^) 
+ 6i(g=')(-4^^, + 7^7,) 

+ b,{q') {g,^ + 2^-^-^q + ^g). (5.5) 

qz 

The above A— current does not satisfy the second Rarita- Schwinger constraint q^J^ = 0. 
This condition would lead to two further relations between the form factors 

(^liq^) = -3a2(g^), (5.6a) 
bi{q') = h{q'). (5.6b) 

17 



However, as will be shown below, the numerical values for these coefficients do not differ 
much from this relation. Consequently, the admixture of spin 1/2 states to the spectral 
function will be small. 

We take the correlator in the euclidean time direction and consider the following 
invariant functions 

nf = ^IV(nf,5^,), (5.7a) 
nf = ^Tr(nJ,V7o), (5.7b) 
nf = ^Tr{U%S^oS.o): (5.7c) 
nf = ^Tr(nf,V^^o7o)- (5.7d) 



As is section 3, we consider the correlators using the vacuum dominance approxima- 
tion. For the four different structures we find 

.2 ^6 



n^ 



TT 



nfo 
nfo 



(5.8a) 
(5.8b) 
(5.8c) 
(5.8d) 



As usual we represent the spectral function for the delta correlator by a resonance 
contribution and a continuum of states starting at an invariant mass sq. The contribution 
of the A-resonance is determined by the matrix element 



<0\jX\^{q,s)> = Aai 



' 2m/\ 



—iqx 



(5.9) 



where u^aiQ, s) denotes a free Rarita-Schwinger vector spinor normalized according to 

q + niA 



^u^{q,s)u^p{q,s) = 



. (5.10) 



2mA V 3 3 3 uia 

The short distance behavior of the coefficients of the spectral function follows from match- 
ing the leading order singularities in the operator product expansion of the correlator with 
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the r.h.s. of eq. (5.5), giving 

<q') = a,{q') = -^ai(g')- (5.11) 

Note that this result is indeed close to the pure spin 3/2 case 02 = — |ai. Using (5.10) 
and (5.11) the invariant correlation functions introduced above are given by 

mq)\ 



nf(r) = 2m^\\D{m^,T) + ^^^ r dssD{^s,T), (5.12a) 

47r^ J so 

Q /-oo 

n^(r) = -2\lD\m^,T)-,^^j^^dss^D'{^s,T), (5.12b) 
n^(r) = -^D'(mA,r)-^ Trf. (-D'(v/^,r) + ^Z)(v/^,r)), (5.12c) 

Tfl/^T ATT J so \T o / 



^ i^-D'iV~s,T) - -D{Vrs,T) + -DWs,T)j. (5.12d) 



1607r4 

The results for the corelation functions are given in figure 4 and table 3. The most 
important difference as compared to the nucleon correlation function is the fact that 
the strong enhancement over the free correlator which we observed in 11^ and 11^ is not 
present in 11^ and Ilf^. This means that the A correlation function is qualitatively different 
from the nucleon one. Note that this feature is not reproduced by the vacuum dominance 
approximation. In particular, vacuum dominance predicts a stronger enhancement in 11^ 
than in 11^. 

Fitting the parameters in (5.11) to the numerically calculated correlators we find 

mA = 1440 ± 70 MeV, Aa = 0.033 ± 0.005 GeVl (5.13) 

With these parameters, we obtain a very good description of all four correlation func- 
tions. The corresponding delta mass is somewhat bigger than the experimental one. This 
situation is similar to the one encountered in the meson sector 0: the delta contains only 
vector diquarks, and, like the rho meson, it also turns out to be too heavy in the RILM. 

As in the case of the nucleon one can speculate whether we really observe a true bound 
state in the spin-isospin 3/2 channel or whether the data can also be explained in terms of 
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the propagation of three essentially uncorrelated quarks. This question is actually more 
serious in the case of the delta, since in this correlator there are no strongly attractive 
diquark channels present. For that we reason we have compared the measured delta 
correlation functions with the simple constituent model introduced in the last section. 
The results are shown in figure 5. Again, the contribution from three uncorrelated quarks 
clearly underestimates the data. We conclude that although there is less attraction in 
the delta channel as compared to the nucleon, the RILM predicts a bound state in the 
spin-isospin 3/2 channel. 



6. Comparison with QCD sum rules and lattice calculations 

Predictions for baryonic correlators based on QCD sum rules were obtained by several 
groups. We show the latest results by Belyaev and loffe pSf , which are compared to those 



from Chernyak et al. (see and references therein). The results of Chernyak et al. are 
based on somewhat different correlators and fitting procedures. In the nucleon case the 
two predictions agree very well (and therefore we have plotted only one of them), but 
they strongly disagree for the delta. 

Chernyak et al. have also performed a more general calculation of 'baryonic wave func- 
tions' and have concluded that there exists a qualitative difference between the structure of 
octet and decuplet particles. Limited data on hard exclusive reactions involving those par- 
ticles seems to confirm this important conclusion. Our results, based on the RILM, also 
show that the octet and the decuplet baryons have very different correlation functions. In 
particular the coupling constants (proportional to the probability to find all three quarks 
at the same point) turn out to be drastically different! 

Our results are compared to QCD sum rule predictions in Fig. 6. In the nucleon case 
the RILM results agree very well with all sum rule predictions. In the delta channel this 
is not the case: the predictions from Chernyak et al. are strongly favored. However, we 
are unable to explain why we are in agreement with certain sum rule calculations but 
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not with others. The whole topic is very confusing, since (as it was demonstrated above 
in the simplest case of diquarks) the agreement of OPE predictions with RILM results 
can depend on the current used, and for some quantities OPE results based on different 
currents do not provide consistent answers. 

In general, it is hardly surprising that these two approaches do not generally agree 
for many particular observables. In fact, the 'vacuum dominance' hypothesis is based on 
the simplistic idea of a very smooth, nearly homogeneous distribution of quark and gluon 
fields in space time. On the contrary, the RILM is an example of a model with strongly 
inhomogeneous fields, with gauge field concentrated in large fluctuations of a particular 
structure, in this case small size instantons. 

Proceeding now to a comparison with lattice calculations, let us first make some general 
comments. 'Traditional' lattice calculations use plane-to-plane correlation function^, and 
those are known to have the following two problems: 

• The measured ratio mj^/mp remains much closer to the 'naive quark model' value 
3/2 rather than the experimental one 1.2, even for the lightest quarks studied (see 
e.g. 
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• The lattice results for the nucleon-delta splitting appear to be smaller than the 
experimental value mA — = 293 MeV. 

It was speculated in that the nucleon, as measured on the lattice, seems to be 
'too heavy', probably because instanton-induced effects were under- represented in those 
simulations, for whatever reason. 

In Fig. 6 we compare our results with recent lattice data |T^, which for the first time 



looked at point-to-point correlators. Good agreement is observed, our RILM results lit- 
erally coincide with the lattice data within the error bars of the calculations! This agree- 
ment is somewhat surprising, because our fit to the RILM results gives mj^/nip = 1±0.1, 
which is smaller than the experimental value whereas our prediction for the splitting 



*And sometimes even a specially chosen 'smeared' sources in order to improve the projection onto the 
ground states. 
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rriA — = 480 MeV is larger than the experimental result. Thus, deviations from ex- 
perimental data in the RILM results have opposite trends as compared to the deviations 
observed in standard lattice gauge calculations. 

Provided a fit to the Negele et al. data gives resonance parameters similar to the 
ones obtained by us, this particular set of lattice measurements seems to show deviations 
from experimental data that are opposite as compared to other lattice works. One may 
speculate that the two sets of lattice measurements emphasize the role of confinement 
and instanton-induced forces differently, because plane-to-plane correlators deal with large 
interquark distances, while point-to-point correlators are also sensitive to smaller ones. 

A controversial subject remains the question whether octet and decuplet baryons have 
a qualitatively different structure. Even a naive constituent quark model including a 



hyperfine interaction leads to significantly different mean square radii, e.g. r/^/r^ ^ 
1.4. If so, the probability to find all three quarks at the same point should be about one 
order of magnitude larger in the nucleon compared to the delta. However, lattice studies 



30|| have been unable to find this effect: it was concluded that the nucleon and the delta 



have essentially the same radius. This question certainly deserves further study. 

One more comment we would like to make concerns the large finite size effects we have 
observed in our calculations. Going from 64 to 256 instantons, (or from a box volume 
of 64 fm'* to 256 fm^) we have seen significant changes in some correlation functions, the 
nucleon correlators among them. Obviously, lattice measurements (which use significantly 
smaller boxes) should be subject to similar modifications. 

Clearly, point-to-point correlation functions are much more sensitive to the dynamics 
of the interquark interactions at small distances. They are also more stable with respect to 
finite size corrections. Further efforts in this direction (using the TERAFLOP or similar 
next-generation projects) should be encouraged. 
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7. Discussion 



Finally, having completed the series of three papers devoted to the RILM, we would 
like to make a few more general comments. 

Although the RILM is a model that was originally designed only to describe the most 
global properties of the QCD vacuum, the gluon and quark condensates, it turns out to 
work very well. Not only does it reproduce small deviations from asymptotic freedom at 
short distances, as it was originally anticipated, but it correctly describes a large number 
of mesonic and baryonic correlators, some of which vary over several orders of magnitude. 
The results obtained imply that even in the absence of confinement the lowest states 
in most channels are in fact bound states of two or three quarks. If we translate the 
measured correlators into resonance masses and coupling constants the results agree with 
experimental data at the 10-15 percent level. 

Does this imply that by fixing the correct magnitude of two condensates, one can 
reproduce all other vacuum parameters as a more or less direct consequence? Absolutely 
not: our experience with interacting instanton ensembles shows that even inside the in- 
stanton framework it is very easy to get drastically different sets of correlators, which 
are completely incompatible with experiment, even if the quark and gluon condensates 
are kept more or less fixed. Thus, the excellent performance of the RILM is indeed a 
non-trivial observation, which probably means it has a lot of truth in it. 

In any case, this model definitely 'outperforms' such widely used theoretical approaches 
as QCD sum rules or the Nambu-Jona-Lasinio model, and therefore applying it to a 
broader set of questions appears to be justified. Among the most fascinating topics to be 
studied are nuclear forces at small distances (propagation of 6 quarks together) , hadronic 
structure and wave functions, multi-point correlators and many others. 

The central question posed at the beginning of this paper has been answered: the 
RILM describes the hghtest baryons as some bound states. 

Moreover, their spin splitting comes out naturally, without any need for additional 
assumptions or parameters. The splitting between scalar and vector diquarks is essentially 
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due to the instanton- induced attraction present in the scalar channel 0, |l^. Although 
we have not considered hyperons in this work, the correct dependence of the instanton- 



induced interaction on the strange quark mass |jT^ allows us to speculate that these 
splittings will also come out reasonable. 

Furthermore, comparison between the calculated nucleon and delta correlators pro- 
vides another important lesson, namely: the octet and the decuplet haryons have com- 
pletely different wave functions, because even the coupling constants (proportional to the 
probability to find all three quarks at the same point) are drastically different. Note 
that this picture of baryons is significantly different from e.g. a 'skyrmion-based' one. 
According to the skyrmion picture, a nucleon and a delta are essentially the same object, 
just rotating with a different angular (and isotopic) momentum. According to our cal- 
culations, they are quite different already at relatively small distances due to the strong 
instanton-induced attraction, which is present in octet baryons but absent in the decuplet. 



8. Conclusions 

We have performed a numerical analysis of diquark and baryonic correlation functions 
in the framework of the random instanton liquid model (RILM). Our main findings can 
be summarized as follows: 

(i) The quark-quark interaction between light quarks is very much channel-dependent, 
and it qualitatively resembles the one in the quark-antiquark channels (with opposite 
parity). In particular, similar to the situation in the pion and rho channels, the interaction 
in the isospin 1=0 (scalar diquarks) and 1=1 (vector diquarks) channels are completely 
different: the scalar channel has a much stronger attractive interaction than the vector 
one, with the correlator reaching a maximum value which is about an order of magnitude 
larger than the free quark correlator. 

(ii) If diquarks are treated as physical particles, the scalar and vector diquarks are 
found to have masses of 420 and 940 MeV, respectively. The non-relativistic description, 
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related to masses of heavy-light baryons (minus the mass of the heavy quark), have 
produced = 760 MeV and e^g = 900 — 1000 MeV, in reasonable agreement with the 
limited experimental information and QCD sum rule calculations. 

(iii) Calculation of all six nucleon and four delta correlation functions have been per- 
formed and analyzed. The octet (nucleon) and decuplet (delta) correlators are found to 
be qualitatively different, with the main difference being due to attractive interaction for 
scalar 1=0 diquarks mentioned above. 

(iv) The nucleon parameters agree well with the predictions from different versions of 
the QCD sum rules, but in the case of the delta the results from Chernyak et al. (see p5| ) 
are strongly preferred. 

(v) Our results were also found to be in agreement with the first lattice data on point- 



to-point correlation functions |T3[ . Moreover, this agreement is surprisingly good, the two 
calculations practically agree inside error bars. 
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Figure Captions 

1. Diquark correlation functions in tlie RILM, measured as a function of distance. All 
correlators are normalized to the perturbative results. The diquark channels are labeled 
by the dirac matrix F defining the current. The solid lines show a fit using the diquark 
model discussed in section 2 and the dotted curves are the prediction from the vacuum 
dominance model. 

2. Correlation functions for A-type diquarks (a-b) and S-type diquarks (c-d). The 
channels are labeled as in figure 1. The solid lines now correspond to the heavy- light 
parametrization discussed in section 3. The dotted and dashed lines in figure (a) corre- 
spond to Grosin-Yakovlev prediction and their OPE expression, respectively. 

3. Nucleon correlation functions in the RILM. Channels are labeled as in section 4 of 
the text. The solid line shows the 'nucleon pole plus continuum' model, while the dotted 
curve is the vacuum dominance approximation. 

4. Delta Correlation functions in the RILM. Curves are as in figure 3. 

5. Comparison of nucleon and delta correlation functions with simple quark-diquark 
models. Figures (a) and (b) show the correlation functions (nucleon) and 11^ (delta), 
respectively. Data and solid curves are as shown in figs. 3 and 4. The dotted curves show 
the constituent model introduced in section 4 of the text. In figure 5a, the quark-diquark 
model is given by the dashed curve. 

5. Comparison with QCD sum rules and lattice data for the two correlation functions 
(nucleon) and (delta). The RILM resuhs are given by the solid triangles, the 
lattice data from [|T^] by the open squares. QCD sum rules predictions from Belyaev and 
loffe are shown as long dashed lines, results from Chernyak et al. as short dashed 
lines. 
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Table Captions 

1. Numerical results from fitting the diquark correlation functions in the RILM with 
a 'diquark resonance plus continuum' model. The parameters are defined in section 2 of 
the text. 

2. Results from fitting the diquark correlators using the non relativistic parametriza- 
tion related to the masses of heavy light baryons introduced in section 3. The two different 
values for the Eq correspond to the results obtained from the F = 7^ and cr^i^ diquark 
channels. For comparison, we also quote the results from two QCD sum rule calculations. 

3. Nucleon and Delta parameters obtained from performing a global fit to the six 
nucleon and four delta correlation functions defined in section 4 and 5. Results are also 
compared to the experimental values for the masses and QCD sum rule predictions for 
the coupling constants. 
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this work 


other information 


comment 


ms 


420 ± 30 MeV 


234 MeV 


NJL model 






rriAv 


940 ± 20 MeV 


824 MeV 


NJL model 






rriT 


570 ± 20 MeV 






9s 


0.225 ± 0.01 IGeV^ 


0.135 ± 0.025 GeV^ 


QCD sum rules |2| 


9av 


0.244 ± 0.010 GeV^ 






9t 


0.134 ±0.004 GeV^ 







table 1 





this work 


other information 


comment 




760 ± 30 MeV 


700 ± 150 MeV 
780 MeV 


QCD sum rules 
QCD sum rules 


M 
Y\ 


1 
1 




890 ± 30 MeV 
1005 ± 30 MeV 


1100 ±200 MeV 
780 MeV 


QCD sum rules 
QCD sum rules 


M 
f>i 


1 
1 


/aq 


0.052 ± 0.005 GeV^ 


0.02 GeV" 
0.0225 ±0.0045 GeV^ 


QCD sum rules 
QCD sum rules 


N 
P 


1 




0.011 ±0.002 GeV^ 
0.044 ± 0.008 GeV^ 


0.06 GeV 
0.0225 ± 0.0045 GeV^ 


QCD sum rules 
QCD sum rules 


|14 

m 





table 2 
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this work 


other information 


comment 


TUN 


960 ± 30 MeV 


939 MeV 


experiment 


Ai 


0.032 ± 0.001 GeV^ 


0.035 ±0.008 GeV^ 


QCD sum rules [19|| 


A2 


0.080 ± 0.004 GeV^ 






^0 


1920 ± 50 MeV 








1440 ± 70 MeV 


1232 MeV 


experiment 


Aa 


0.033 ±0.005 GeV^ 


0.050 ± 0.013 GeV^ 


QCD sum rules [IS] 






0.035 GeV^ 


QCD sum rules 


Eo 


1962 ± 101 MeV 







table 3 
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